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ABSTRACT 

Observations of relativistic jets from black holes systems suggest that particle acceleration often 
occurs at fixed locations within the flow. These sites could be associated with critical points that 
allow the formation of standing shock regions, such as the magnetosonic mo dified fast point. U sing 
the self-similar formulation of special relativistic magnctohydrodynamics by iVlahakis fc Konigll we 
derive a new class of flow solutions that are both relativistic and cross the modified fast point at a 
finite height. Our solutions span a range of Lorentz factors up to at least 10, appropriate for most 
jets in X-ray binaries and active galactic nuclei, and a range in injected particle internal energy. A 
broad range of solutions exists, which will allow the eventual matching of these scale-free models to 
physical boundary conditions in the analysis of observed sources. 

Subject headings: acceleration of particles — ISM: jets and outflows — methods: analytical — MHD 
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1. INTRODUCTION 

Jets have been observed around a large variety of astro- 
physical objects, such as young stellar objects (YSOs), 
accreting white dwarfs, X-ray binaries (XRBs), and ac- 
tive galactic nuclei ( AGN) , and are also thought to drive 
gamma-ray bursts (GRBs). In YSOs jets facilitate an- 
gular momentum transport, allowing the central star to 
accrete more matter, and likely play a similar role in 
accreting black hole systems, where in AGN they are 
thou ght to also affect the evolution of their host galaxy 
(e.g. iBest et"aL|[2006t) . While there are enormous dif- 
ferences between the types and scales of objects around 
which jets can occur, the origins of jets seem to be re- 
markably similar, requiring the basic ingredients of in- 
falling/collapsing, rotating matter, and magnetic fields. 

Despite these seemingly simple initial conditions, there 
are currently many outstanding problems in our under- 
standing of jets, from their creation to their matter con- 
tent and internal physics. One important facet of jets 
observationally is their hallmark synchrotron emission 
dominating the radio bands in particular, extending up 
to at least the near infrared (NIR) in XRBs (jFender et al.1 
11997ft . At higher frequencies the typical power-law spec- 
tral energy distributions (SEDs) are generally interpreted 
as optically thin synchr otron emission from acc elerated 
particl es, both in A GN (jMarscher fe Gearl[l985ft and in 
XRBs ijFenderl [2002ft . An effective way to accelerate ra- 
diating particles into a power -law distribution is via dif- 
fusive shock acceleration (e.g. lBeUll978t lDrurvill983ft off 
scattering centers in turbulent plasma flows. Once initi- 
ated, this process must to be distributed throughout the 
flow to account for the lack of spectral aging over vast 
distances along the jets (e.g. lJester et al.ll2001| ). 

But where does the acceleration itself begin, and what 
triggers it? There is increasing evidence that the start 
of this region is offset from the central compact object. 
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For example in the jet of the AGN M87, recent VLBI 
observations show the synchrotron emi ssion starting at a 
region offset by ~ 1 00r n from the core (|Junor et alj|1999t 
IWalker et al.ll200"8ft . Similarly, the start of power-law ac- 
celeration in compact jets would be indicated by a tran- 
sition from optically thick emission, with a flat/inverted 
spectral index, to an optically thin power-law at a dis- 
tinct location in the SED. Such a break has been observed 
directly so far only in one source, the Galac tic XRB GX 
339-4, in the NIR ([Corbel fc Fended 12001 . during the 
"hard" or "nonthermally dominated" accretion state as- 
sociated wiHi_cojirr>aj3tje^ state definitions 
in, e.g.. iMcClintock fc Remillardll200"6ft . Because of the 
stratification of emission regions in compact jets (e.g. 
iBlandford fc Koniglll979D . the lower the frequency where 
this turnover occurs, the further the location along the 
jet where particle acceleration starts. In XRBs, a break 
in the NIR corresponds to an offset of ~ 10 — 1000 r g , 
and models of the broadband data of most black hole 
XRBs in t he hard state so far seem to require such a 
break (e.g. iMarkoff et al.ll200lLl200aiNowak et al.M 2005t 
iMigliari et alJl2007b iMaitra et alJl2009ft . 

The fact that the start of the acceleration region seems 
to occur at roughly the same location in several sys- 
tems could be indicative of a critical point occurring in 
a magnctohydrodynamical (MHD) flow, particularly the 
magnetosonic modified fast point (MFP). At the MFP 
the collimating magnetic field lines turn inwards towards 
the jet axis, potentially leading to rccollimation shocks, 
while at the same time the flow becomes causally discon- 
nected so shocks can occur without disrupting the flow 
upstream. Such a shock region thus would occur at a 
fixed location in the flow, closely connected to the MFP, 
and would be an ideal location for particle acceleration 
to begin. We wish to investigate the feasiblity of this 
premiss in this paper. 

Because of the complexity involved in relativistic MHD 
including strong gravity, many groups are using the re- 
sults of simulations to study the formation and develop- 
ment of jets. These simulations often show the develop- 
ment of a steady outflow in which the magnetic field is re- 
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markably self-similar and axisymmctri c near the launc h 
point (see, e.g., Fig. 2 and Fig. 11 in iMcKinnevj 120061 ). 
However to study specifically the development of criti- 
cal points in the flow and their dependence on external 
boundary conditions, current MHD simulations either do 
not extend far enough from the black hole, or if they 
do, they are too computationally expensive. Assuming 
that the resulting flows retain a self-similar structure, 
at least when gravity does not dominate as indicated in 
the simulations, w e adop t the formalism developed by 
IVlahakis fc Konigll ()2003L hereafter VK03). By assum- 
ing axisymmetry and a self-similar field line geometry, 
VK03 reduce the exact equations of special rclativistic 
MHD to a onc-dimensional problem. Although VK03 fo- 
cused on jets in GRBs, this treatment is also applicable 
to other MHD jets such as in AGN and microquasars. 

In an ea rlier self-simila r trea tment that was non- 
rclativistic, IVlahakis et al.1 (|2000l ) presented a solution 
where the flow crosses the MFP at a finite height above 
the disk. In VK03, however, they only found relativis- 
tic solutions with an MFP occurring at infinity (meaning 
that the flow asymptotically approaches a perfect cylin- 
drical geometry). In this paper we extend the study of 
VK03 and derive new solutions in which the relativistic 
flow crosses the MFP at a finite location above the disk. 
In § 2 we describe the VK03 model and our method for 
exploring the full parameter space of solutions. In § 3 we 
present the first relativistic solutions that pass through 
an MFP. In § 4 we discuss our results, and the depen- 
dence of the MFP location on the model parameters. Wc 
also describe how this work sets the stage for further de- 
velopment to connect the flow to regions near the disk 
where gravity can no longer be ignored. § 5 contains our 
conclusions. 

2. METHOD 

2.1. Background 

During the acceleration and collimation of jets, 
magnetic fields are thought to efficiently extract 

rotational energy fr om either the compact object 

Blandford fc Znaiekl Il977f ) or the accretion disk 



Blandford fc PavndllQSir The latter models are in the 



Newtonian limit, with the matter considered cold, mean- 
ing there is negligible thermal pressure causing bulk ac- 
celeration to non-relativistic velocities. 

These cold, non-relativist ic solutions we re generalized 
to the relativistic regime by iLi et al.l (|1992[ ) , allowing the 
bulk velocity to attain relativistic speeds. VK03 further 
extended the solutions to include the "hot" regime, al- 
lowing the random motions of the particles to become 
relativistic and thus the jets to be hydrodynamically ac- 
celerated even at the base where temperatures are high. 
It is this last scenario that we base this work upon. 

Starting from the equations of time-dependent special 
relativistic MHD, we make the following assumptions to 
render them more tractable: ideal MHD, no gravitational 
field or external force (and thus self-similarity holds), 
axisymmetry, a zero azimuthal electric field and time in- 
dependence. Following the terminology of VK03, after 
scaling the equations to make them non-dimensional, we 
are left with two coupled differential equations (equations 
(|A4[) and (|A5[) in the appendix). Combining these two 
coupled differential equations, we obtain a single equa- 



tion for dM 2 /d# (equation (TA"9|) . with M the Alfvenic 
Mach number and 9 the angle of the point on the field 
line with the axis of symmetry) which acts like a "wind 
equation" , much ak in to the win d equation of the Parker 
solar wind model (Parker 1958). This wind equation, 
along with the other algebraic equations (see appendix) , 
can be solved for the velocity, magnetic and electric field 
strength, density and pressure along a field line. Due to 
the self-similar assumption, once solutions are obtained 
for one field line, all other field lines can be obtained by 
simple scaling. An example of this self-similarity and the 
meaning of some of the parameters used can be found in 
Figure [TJ 

Instead of the single critical (or sonic) point of the 
Parker solar wind model, due to the inclusion of magnetic 
fields, the obtained wind equation has three locations 
where the denominator crosses zero. Starting from the 
accretion disk, these are the modified slow point (MSP), 
the Alfven point, and the MFP. The MSP and MFP 
are also called the slow and fast magnetosonic separatrix 
surfaces. The Alfven point is the location where the rel- 
ativistic collimation speed of the flow towards the axis 
(Vg) is given by 



(jv e y 



Bl 1 



(1) 



where 7 = 1/(1 — V^/c 2 ) 1 / 2 is the Lorentz factor, B the 
strength of the magnetic field, x the cylindrical radius 
in units of the light cylinder radius, po the baryon rest- 
mass density, and £c 2 the specific relativistic enthalpy 
(the variables are described in more detail in H2.2\i . The 
denominator of dM 2 /d# with the Alfven point divided 
out, can be expressed as 



c 

C?Bf(l 



1V9 
c 

™2\ 



El 

„2 



B 2 -E 2 



47T/?o£c 2 



with 



c 2 iirpo^c 2 



V 2_„2 (r-m-i) 



(2-r)f + r-i' 



(2) 
(3) 



c the velocity of light, E the strength of the electric field, 
and r the polytropic index. The MSP and MFP are, by 
definition, the locations where & = (VK03). 

At every critical point, the numerator of dM 2 /d6> 
should also pass through zero to ensure a smooth cross- 
ing. This translates into a regularity condition at the 
critical points and fixes the value of a free parameter. 
Even though the MSP should be crossed smoothly to 
obtain a solution that describes the entire jet from the 
accretion disk to the termination point, gravitational ef- 
fects cannot be ignored at the MSP. Since the equations 
do not include gravity (as it is not compatible with the 
self-similarity assumption in relativistic flow) , we do not 
try to fit for the MSP. Therefore we fit two critical points, 
and correspondingly two parameters are fixed, in our ap- 
proach described below pj± for the Alfven point and ctm 
for the MFP. 

The physical importance of the MFP is that it is the 
location where not even the fastest signals can travel up- 
stream anymore, meaning anything downstream from the 
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Figure 1. Sketch of self-similar field lines projected onto the 
meridional plane. For any 8 the values of the variables describing 
the field line are exactly the same, only scaled by their respective 
distances. 

MFP is causally disconnected from the region upstream. 
If, for example, a shock were to exist beyond the MFP, it 
could not disrupt the flow leading to that shock, allow- 
ing it to be a permanent feature. As mentioned above, 
at the MFP there is a component of the velocity heading 
towards the axis. This can lead to a collimation shock 
shortly beyond the MFP, causing the magnetic energy 
to be converted to particle energy and the jet to become 
kinetically dominated. Another possibility is that the 
flow remains magnetically dominated, and, after reaching 
a minimum radius, bou nces back, retaining an ordere d 
magnetic field (see, e.g.. iContopoulos fc Lovelace] 1 19941 ). 
This is not in conflict with the statement in VK03 that 
the only physically acceptable case in the super- Alfvenic 
regime is for the flow to become asymptotically cylin- 
drical, as this statement only applies to solutions with 
F > 1. 

2.2. Model parameters 

The prescription we are following from VK03 has 9 free 
parameters that determine the solution, whose effects are 
described below. Following the same notation, a Roman 
subscript A signifies the value of a variable at the Alfvcn 
point and an Italic A denotes a value with respect to the 
poloidal magnetic flux function. 

2.2.1. Free parameters 

The exponent F determines the current distribution. 
A value F > 1 corresponds to the current-carrying 
regime, with higher values of F ensuring faster collima- 
tion, but if F < 1 we are in the return-current regime. 
The restriction on F is that it cannot be negative: F > 0. 
Although we consider F to be a free parameters, for this 
paper we chose to keep it fixed at 0.75. 

The adiabatic index T can have values of 4/3 for rela- 
tivists and 5/3 for non-relativistic solutions. 

9 A gives the angle where the Alfven point is located 
with respect to the axis of symmetry. It is limited by the 
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Figure 2. 3D plot of solutions with an MFP with parameters 
F = 0.75, T = 5/3, x\ = 0.75 and 6» A = 50°. On the front side 
the solid blue and dasned red surfaces are connected. The lower 
surface does not extend all the way to the -i/>A-axis for all values 
of ip A . Because we favor solutions with strong magnetic fields, we 
focus on the upper surface, with high ctm- 

(A color version of this figure is available in the online journal.) 

value of ip A : 9e (90° - ip A , 90°]. 

i/j a gives the poloidal slope of the field line at the 
Alfvcn point with respect to the accretion disk. This 
in turn is limited by 9 A : i>A e (90° - 6 A , 90°]. 

x\ is the radius squared of the Alfven point in terms 
of the light cylinder radius. For x A — >• 1 the solu- 
tion becomes more force-free. The allowed values are 

4g (0,1). 

<tm is th e magn e tizati on parameter in the monopole 
solution of iMichel (|1969l) and is related to the mass-to- 
magnetic flux ratio. The constraint is ctm > 0. 

q is the dimensionless adiabatic coefficient and is con- 
stant along a field line. For a large value of q, the specific 
rclativistic enthalpy of the matter is high, for q — > equa- 
tion (|A2|) shows £ — > 1 and the flow is cold. Therefore 
q > 0. 

p A is the derivative of M 2 with respect to the polar 
angle 9 at the Alfven point. Accelerating flow implies 
p A < 0. 

Bqvj 2 ~ f , with reference magnetic field Bo and refer- 
ence length wo, is the scaling of the solution, relating the 
dimensionless values to physical dimensions. We do not 
yet apply our solutions to specific black hole systems, so 
this parameter is not used here but it will be important 
for future applications of our solutions. 

The smooth crossing of the Alfvcn point is ensured by 
calculating p A from the corresponding regularity condi- 
tion, given by equation (|A8[) . In the same way we de- 
termine cm by crossing the MFP, although this param- 
eter sometimes can have two values (see Figure [2]). This 
leaves F, T, 9a, tJja, x\, and q (and Boz&q~ f ) to satisfy 
the boundary conditions at the source and, indirectly, at 
the end of the jet. 

2.2.2. Other parameters and variables 

There are also parameters derived from the above val- 
ues: 

/i determines the total energy-to- mass- flux ratio (/ic 2 ) 
and is conserved along a field line. This parameter is 
determined from equation (|A7[) and /j, > 1. 
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£7 A is the value of the magnetization function ct, de- 
fined as the Poynting-to-matter energy flux, at the A lfve n 
point. This parameter is determined from equation (|A6|) . 
As /i cannot be negative, from equation (|A7|) follows 

And finally wc describe the other variables used in the 
equations: 

£ determines the specific (per baryon mass) relativistic 
enthalpy (£c 2 ). If £ = 1 we have cold, pressureless mat- 
ter. £ will drop from the high temperatures at the base of 
the jet as matter is mainly accelerated hydrodynamically 
and at some point above the disk drop down to 1. From 
this point on all acceleration is magnetic. This variable 
is determined from equation (|A2[) . 

M is the Alfvenic Mach number, the velocity of the 
flow in terms of the Alfven velocity. 

G is the radius in terms of xa and is therefore equal 
to 1 at the Alfven point. 

2.3. Numerical method 

To find solutions with an MFP we obtained expressions 
for dM 2 /d9, given by equation (|A9[) . and for dtp/dO by 
combining the derivative of the Bernoulli equation with 
the transfield equation using the determinant method. 
Because dip/d6 is very unstable near the Alfven point, 
as the numerator has a first order zero point there and 
the denominator a second order one, we reverted to the 
Bernoulli equation, given by (|A3[) . to determine ip. To 
start off the integration from the Alfven point we specify 
F, r, 6a, x a , "0a, q and an initial guess for ctm. We deter- 
mine pa from the Alfven regularity condition, equation 
(|A8|) . Integrating outward from the Alfven point, we de- 
termine whether the numerator or denominator crosses 
zero first and adjust ctm accordingly until both cross at 
the same time. We then proceed to explore the range 
of solutions which cross both the Alfven point and MFP 
(see Figure EJ). For plotting purposes, we have divided 
out the factor a; 6 (l — M 2 — x 2 ) 2 in both the numerator 
and denominator. 

As discussed above, because there is no gravity in the 
model, we do not try explicitly to cross the modified slow 
point (MSP). Gravitational effects should play a large 
role close to the black hole, and the self-similar equations 
cannot predict accurately where the MSP is located. 

2.4. Approach 

To begin our exploration of parame ter space, we chose 
a solution from Vlahak is et al.1 (|2000( ) known to have an 
MFP (specified below their Figure 4) with parameters 
x = 0.75, 7 = 5/3, 0» = 60°, V* = 45°, k 2 = 15, 
H = 10.9239, A 2 = 2.79 35, p« = -5-5744 and e = 9.4487. 
By comparing terms in I Vlahakis et al.l ()2000l ) and VK03, 
it is possible to translate these parameters to the param- 
eters used in VK03. Because we are using the relativistic 
equations from VK03, the parameters of our first solu- 
tion (see Table [1]) differ slightly from the corresponding 
parameters above and we vary ctm to obtain a critical so- 
lution again. From this solution we were able to traverse 
parameter space, while allowing only critical solutions 
with an MFP. By increasing x\ we were able to obtain 
higher velocities of the jet. After achieving relativistic 
velocities, we focused on finding solutions with higher 
values of q, as the solutions so far were cold. But for 
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Figure 3. Solution c. For parameters see Tablefl] The scaling for 
the y axis is a combination of a linear and logarithmic part, using 
the function sign(x) log 10 [l+abs(x)/10 — 12 ] . The vertical black line 
gives the location of the Alfven point (at 8 « 0.87). Even though 
the numerator (solid red line) and denominator (dashed green line) 
change sign, it can be seen from their ratio (short dashed blue 
line, dM 2 /dd) that the MFP is crossed smoothly. We ceased the 
integration shortly after the MFP. 

(A color version of this figure is available in the online journal.) 

a fixed value of x\ there is a maximum value of q that 
produces critical solutions crossing the MFP. This is due 
to the fact that the collection of all solutions form a sur- 
face in the multidimensional parameter space, and we 
had reached a maximum for q for the fixed values of the 
other parameters (see Figure [2]). 

The surfaces of valid solutions have roughly the same 
appearance for the explored range of 9a and x\- To 
describe the effect of 9a and x\ on the solutions, we 
approximate the graph as a cone with the base in the 
%pA, OM-plane and with the maximally allowed value of 
q as its height. If we increase a; a the base of the cone 
shrinks while moving to higher ctm, and the height de- 
creases. The latter limits the maximum value that x\ 
can be increased to. If we decrease 9a, the height in- 
creases, indirectly allowing higher values for x A . The 
area of the base becomes bigger and shifts to higher ipA 
and ctm, with the upper surface becoming steeper. Due 
to the shape of the surface, it is possible to increase any 
two parameters of x\, q and ctm at the expense of the 
third. 

By extending our search to three parameters ("0a, <? 
and ctm), we were able to move around this point and 
continue increasing q. This revealed a multidimensional 
surface that is double- valued in ctm- An example of this 
surface is shown in Figure [H which also includes our most 
relativistic solution presented below, solution c. The nu- 
merator, denominator, and acceleration of M in this 
latter solution are also plotted in Figure [U 

3. RESULTS 

In this section wc present the various solutions crossing 
the MFP, from the first one that we found to one with 
relativistic temperature and bulk flow we sought, while 
describing the features particular to a certain solution. 
The parameters of our solutions are given in Table [T] 
and the main properties in Figure [H As we are not yet 
applying our solutions to specific black hole systems, we 
do not use the scaling parameters BqWq~ f and wq. 
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Table 1 

Parameters of solution 



Solution 


F 


r 


T 2 
A 


q 


8a 
(deg) 


V'A 
(deg) 


°"M 


PA 


fA 




a 


0.75 


5/3 


4.7676 x 10 -3 


2.5348 x 10" b 


60 


45 


2.14228 X 10" 4 


-5.54314 


5.42665 x lO -3 


1.01197 


b 


0.75 


5/3 


0.75 


2.5 x lO" 6 


60 


45 


1.01241 


-1.50363 


3.14708 


5.44084 


c 


0.75 


5/3 


0.75 


0.12 


50 


55 


2.53981 


-1.66651 


2.6039 


9.85117 



Note. — The values for the first six parameters (F through ?/>a) are exact, for the last four (ctm through /i) they are rounded off. 



Solution a 



Solution b 



Solution c 




10 100 1000 

G 



Figure 4. Properties of solutions. Row 1 shows the geometry of the field line, where the height has the same scaling as the cylindrical 
radius. Row 2 shows the Poynting-to-mass flux ratio (S = — toJIB^/^ac 2 ), the Lorentz factor 7, and the specific relativistic enthalpy. 
Row 3 shows the causal connection opening angle arcsin(l/7), and the opening half-angle of the outflow (it/2 — ip), which is negative after 
overcollimation. Row 4 shows the squares of the Alfvenic Mach number (M 2 ), and the radius in units of the light cylinder radius (oc 2 ). 
(A color version of this figure is available in the online journal.) 



3.1. Solution a: A cold, slow jet 

This solution is the c losest to the non-rela tivistic pa- 
rameter values given in IVlahakis et all (j2000T ) that con- 
forms to the Alfven regularity condition (ARC), which 
is the transfield equation at the Alfven point. The so- 
lution crosses the MFP at 6 re 0.15 or 8.6° (G re 15.3). 
The Alfven point is located at 6 re 1.05 or 60°. The top 
left panel of Figure [4] gives the meridional projection of 
the magnetic field lines. The jet overcollimates after a 
maximum radius of almost 16 times the Alfven radius at 
6 re 0.20 or 11.5°, shortly before the MFP. The second 
left panel of Figure [4] shows the flow is cold throughout 



(£— 1 <C 1) due to the very small value for q and low x\. 
As M 2 



1 — x A at the Alfven point, equation (|A2j) gives 
a value for £a very close to 1. The energy of the mat- 
ter £ (including the dominant rest mass energy) is much 
higher than the energy in the magnetic field throughout. 
Therefore, even though the magnetic acceleration is ef- 
ficient, the jet is not accelerated to relativistic velocities 
(7 < 1.02). The third center panel shows the "causal 
connection" opening angle arcsin(l/7) and the opening 
half-angle of the outflow, which goes from 60° to a few 
degrees overcollimation. Although the causal connection 
opening angle has little importance for non-relativistic 
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flows, as it remains very close to 90° since 7 ~ 1, it is 
shown for completeness. 

3.2. Solution b: A cold, fast jet 

After our first solution, we increased the velocity of our 
jet by increasing x\. As the top center panel in Figure 
[4] shows, after a long period where the field line remains 
almost parabolic, the jet in this solution overcollimates as 
well. This is caused by magnetic hoop stresses and may 
allow a shock region to develop beyond the MFP. The 
Alfvcn point is again located at 9 ~ 1.05 and the MFP 
at 6 « 0.063 or 3.6° (G » 145). The second center panel 
shows the Lorentz factor and the enthalpy of the flow. 
The flow here also is seen to be cold (£ ~ 1), meaning 
the jet is mainly magnetically accelerated, which is again 
due to the small value of q. The Poynting-to-mass flux 
ratio (S = —tijVLB^/^ac 2 ) decreases, showing magnetic 
energy being transferred into kinetic energy, with the 
flow reaching a Lorentz factor of 2.8. The bottom center 
panel shows the squares of the light cylinder radius, x = 
wVt/c, and the Alfvenic Mach number, M, When x 2 = 1 
the light surface is reached, which is the radius where the 
field circular velocity reaches the speed of light. 

3.3. Solution c: A warm, very fast jet 

After having achieved a relativistic solution for the cold 
plasma case, we would like to find solutions with an in- 
creased flow temperature. To do so requires increasing 
the value of £. As M 2 is given by 1 — x\ at the Alfvcn 
point, equation (|A2|) shows that increasing x\ and/or q 
has the desired effect. Unfortunately, for larger x\ the 
maximum value of q decreases. By choosing a lower value 
for 9a the attainable values for x\ and q are increased, 
leading to a warmer flow. This solution is shown in the 
third column of Figure |4] The Alfven point is located 
at 9 « 0.87 or 50° and the MFP at 6 m 0.041 or 2.4° 
{G « 87.4). Near the beginning of the flow £ » 2.9. 

The Lorentz factor of the flow at the MFP is 8.3, which 
is mainly due to the high initial Poynting flux. 

It can be seen in the second right panel of Figure[4]that 
£ always drops to 1 and from the fourth right panel that 
M 2 always dominates x 2 near the MFP. As the Lorentz 
factor is given by 

tl l-M 2 -xj 

£ 1 - M 2 — x 2 V 1 

the final Lorentz factor is approximately ji. This means 
that while the jets may start as Poynting flux-dominated, 
eventually they convert most of their Poynting flux and 
become kinetic energy-dominated. 

3.4. Location of the MFP 

Since we are interested in the location of the MFP, we 
would like to know how it depends on the model param- 
eters. As it is challenging to sample the full parameter 
space, we will focus on the region around the solution 
closest to observed jets, solution c. By allowing all pa- 
rameters to vary and looking at the effect this has on 
the location where the MFP occurs, we can draw the fol- 
lowing conclusions: the MFP moves outward (smaller 9) 
when the Alfven point occurs at a smaller angle, (lower 
9a), when the temperature at the base of the flow is in- 
creased (higher q), when the flow at the Alfven point 



is already close to collimation (higher iJja), or when the 
Alfvcn point moves closer to the light cylinder radius, 
making the flow more force- free (higher ia)- 

4. DISCUSSION 

We have succeeded in obtaining new, solutions for a 
relativistic, magnetized flow that smoothly crosses the 
MFP at a finite height (9 > 0) above the system equa- 
tor. These solutions suggest that it should be possible to 
construct better, more MHD-consistent jet models where 
the location of the acceleration region is determined a 
priori from the physical boundary conditions. 

So far none of the solutions derived remains Poynting 
flux-dominated up to the MFP, which is probably due to 
the relatively small value of cm found so far. Increasing 
x\ and especially decreasing 9a will allow higher values 
of <tm to be used. The same can also be done at the 
expense of q. 

We also have described some of the relations between 
the different parameters (see §2.2j) and the effect they 
have on each other. Having a more force-free solution 
(higher x A ) decreases the allowed range of temperatures 
(q) and collimation at the Alfv'en point (V'a), but at 
the same time allows a higher value of cm that provides 
a critical solution. Keeping all other parameters fixed, 
there is a maximum value of x\ for which a solution is 
possible at all. This maximum may be increased by mov- 
ing the Alfvcn point closer to the disk (smaller 9a)- This 
change has the additional effects of allowing a broader 
range for the collimation angle at the Alfven point (^a) 
while at the same time shifting this range towards higher 
collimation. It also allows for a higher temperature of the 
flow (q) or a higher magnetic field strengths (ctm)- Any 
pair of parameters x\ , q and ctm may be increased at the 
expense of the third. 

Two parameters not varied so far are F and T. Higher 
values for F should ensure faster collimation and might 
therefore be very important for the exact location of the 
MFP. Similarly, we may want to explore the T = 4/3 case 
for jets with extremely relativistic temperature. How- 
ever, for the weaker jets in AGN and XRBs that we plan 
to target, the radiating particle distributions are gener- 
ally thought to peak at mildly relativistic energies. 

5. CONCLUSION 

If the start of the particle acceleration region in steady 
jets is indeed associated with the magnctosonic fast crit- 
ical point in the bulk flow, then our results support the 
conclusion that such a region could occur at a fairly sta- 
ble location about the launch point. All of the solutions 
found overcollimate shortly before the MFP, as would 
be expected for the initiation of shock development. By 
starting with a sol ution that crosse s the M FP in the non- 
relativistic case of lVlahakis eTall (f2000T) . we were able 
to extend the solution through the multidimensional pa- 
rameter space towards relativistic velocities and temper- 
atures, while retaining the critical point. This feature 
sets our results apart from the work of VK03, whose for- 
mulism we adopted, in that they are immediately appli- 
cable to observed compact jet sources with an optically 
thick-to-thin break in the synchrotron spectrum, such as 
hard state XRBs and weakly accreting AGN. Our most 
promising solution is a jet outflow with mildly relativistic 
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temperature, and Lorentz factor of ~ 10, also appropri- 
ate for the steady jets in both XRBs as well as AGN. 

It is clear that a wide range of parameter space is 
left still unexplored, that can be exploited for matching 
physical boundary conditions appropriate to known as- 
trophysical sources. However, before a radiative model 
can be constructed around the dynamical "backbone" 
provided by the solutions presented here, a prescription 
for including gravity must be included, to extend these 
solutions through the MSP and allow connection with 
a physical model of the accretion flow/corona. We are 
currently working on matching these necessarily non-self 
similar solutions to those presented here, which will be 
presented in a separate work. Once this solution is in 
place, we will have a mu ch more physically consistent 
model (compared to, e.g., [Markoff et al. 2005) to use in 
the fitting of data from accreting black holes across the 
mass scale, which show compact, steady jets. As there 



seems to be no shortage of possible solutions, we are con- 
fident that we can match physical boundary conditions 
with critical solutions. 
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APPENDIX 



EQUATIONS 

Here we list the equations we have used for reference. See tj2.2l for a description of the parameters and variables. 



x = xaG 



(Al) 



Al = q- 



(£_l)i/(r-i) 
The Bernoulli equation is given by 

I? G 4 (l - Al 2 - x\) 2 - x 2 (G 2 - Al 2 - x 2 ) 2 
£ 2 G 4 (l - M 2 - x 2 ) 2 

dG 2 2G 2 cos(V) 



F 2 a 2 A M i sm\6) 
£ 2 x 4 cos 2 W> + 6>) 



The transfield equation is given by 

d 



Gsin 2 (<9) 



df) 



tan(V-' + 6) 



d(9 



1-M 2 



sin(6>) cos(i/' + 0) 



(A2) 

(A3) 
(A4) 



G 



l-G 2 



-sm z (9) 



M 2 + Fx 2 - F + 1 



X p^^JL X 



cos 2 (V> + 9) 



G - M 



„2\ 2 



T-lF-2i{i~l)x A 



T F 2 c 



M 



M 2 



The magnetization function a and the fractions at the Alfven point are given by 

1x\ cos(-0a) / 1 — M 2 — x\ \ 



a A = 



PA sin(0A) cos(Sa+i/'a) ' 



( 1-G 2 \ _ ta/x\ 



(- 



M 2 -x 2 



-Af-xi 



CTA + 1 ' 

2: a (<ta + 1) ' 



(A5) 



(A6) 



We obtain fj 2 by inserting these relations into the Bernoulli equation. 



A' 2 = 



(cta + I) 2 



<T A (1-X A )Y 



2 t-2 



F 2 a 2 M (l-x 2 A y S in 2 (6 A ) 



x\ cos 2 (6» A + -0a) 



(A7) 



The Alfven regularity condition is obtained by substituting these relations into the transfield equation. 



F 2 <(l-xl)(a A + l)"sin( 
cos 2 (6> A +i>A) 



- 2 



r-i(F -2)(q-i)(i- a; |) 



sin(^A) + 2 cos^a) sin(0A + ?/>a) 



0"A 



+ ^ [(F-l)(l-,i)-l] 



= [x 2 A - a A (l ~ x 2 A )] 2 
The wind equation is given by 



(F- 1)4(1 -4) -2 



r - i 



where 



AM 2 dB 2 - C 2 B l 



d(9 



A 1 B 2 -A 2 B 1 



(f-2)^{4-[4-- A (i-4)] 2 } (A8) 

(A9a) 



Ax 



+ G 2 {1- M 2 - x 2 ) 



^(r+i) g r-i 



, F 2 <M 2 sin 2 (g) 

x 4 cos 2 (V' + 6<) ^2r c^r-2^j 



Si 



-2 J F 2 cr 2 4 M 4 sin 2 (e) tan(?A + 9)G 2 (1 - il/ 2 - x 2 ) 2 
,dG 2 



Ci = x 4 cos 2 {i> + 6)n 
AG 2 



dO 



(1 - M 2 - xlY + 2x 2 A (G 2 - M - x 2 ) (1 - 4) 



+ -r^-(l - M 2 - x 2 ) (1 - M 2 - 3a; 2 ) [£ 2 x 4 cos 2 (V> + 9) + F 2 al l M 4 sin 2 



(16* 

+ 2F 2 4A/ 4 sin 2 ((9)G 2 (l - M 2 - x 2 ) 2 
A 2 = - sir 2, 
B 2 = shi 2 (0) 



cos(6>) 1 dG 2 



sin(6>) G 2 d<2 



tan(i/' + #) 



in 2 (0) tan(^ + 0) 
1 - M 2 - x 2 



cos 2 (-0 + 0) 



G 2 = -sin 2 (0) 



1 - M 2 



cos 2 + 6) 
+ {F-l) 

4 2 
X fJi X 



sin 2 (6')tan('!/' + 9) 

2 



1-G 2 



4 2 2 

x A M a; 

^2-2 

^2 / G 2_ il/ 2_ X 2N 2 



sin 2 (6>) 



2 dG 2 , , r2 2 , 1 dG 



ns M 2 + Fx 2 — F +1 



F 2 a 2 A M 2 \ 1 - M 2 - x 2 



cos 2 (^ + 9) 
T - 1 F - 2 £(£ - l)x 4 



T F 2 c 



M 



M 2 



(A9b) 
(A9c) 



(A9d) 
(A9e) 
(A9f) 



(A9g) 
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